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A GENERALIZED REPRESENTATION FORMULA FOR 
SYSTEMS OF TENSOR WAVE EQUATIONS 

ARICK SHAO 



^_ ■ Abstract. In this paper, we generalize the Kirchhoff-Sobolev parametrix of 

J^ ' Klainerman and Rodnianski in 1121 to systems of tensor wave equations with 

additional first-order terms. We also present a different derivation, which 
better highlights that such representation formulas are supported entirely on 
past null cones. This generalization of [12] is a key component for extending 
Klainerman and Rodnianski's breakdown criterion result for Einstein-vacuum 
spacetimes in 1141 to Einstein-Maxwell and Einstein- Yang-Mills spacctimcs. 



< 

^ . 1. Introduction 

-i— > . 

Let (M, g) denote a (1 + 3)-dimensional time-oriented Lorentzian manifold, with 
Levi-Civita connection D and Riemann curvature R. Also, fix an integer n > 0, 
along with n integers r (1) , . . . r ( ™> > 0. For each m,c £ {1, . . . , n}, we let 3> <m) and 
ij <m) denote tensor fields on M of rank r <m) , and we let P< mc > denote a tensor field 
£> ■ on M of rank 1 + r (m) + r (c) . u Assume that these objects satisfy the system 

ov 
^. : 

\f^ . of tensor wave equations on M, where D g = g a ^D a p is the covariant wave operator, 

and where / and J are collections of r <m) and r <c) indices, respectively. □ Note that 
if n — 1 and P (11) vanishes, then ([1} reduces to the tensorial wave equation 

(2) D 5 $ = *. 

The aim of this paper is to generalize the representation formula of |12j , which 
handled the setting ([2]), to also treat systems of the form (JTJ). In other words, we 
determine a formula for the $ <m) 's at a point p € M in terms of the <I> (m) 's and 
, 3> <m) 's along a portion of the regular past null cone about p. In addition, we present 
a derivation different in nature from that of [12] . which allows us to weaken the 
assumptions required for this formula to be valid. Lastly, like in [12] , we will extend 
this representation formula to arbitrary vector bundles. 

1.1. Prior Results. The model problem for (flj is the scalar wave equation 

(3) D0 = V, </>,V>eC c 
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2 ARICK SHAO 

on the Minkowski spacetime R 1+3 . Suppose 4> solves ((3]), with initial data 

cb\ t=0 = a e C°° (M 3 ) , 9 t <£| t=0 = ai er (R 3 ) . 

From standard theory, cf. Sec. 2.4], at a point (t, x) £ (0, oo) x M 3 , we can 
express 4>{t, x) explicitly in terms of the initial data «o and o.\. 

Recall that we can decompose <f> = <f>\ + fa, where fa satisfies Q but with zero 
initial data, while fa satisfies the homogeneous equation Dfa = with initial data 
a® and a±. The term fa can be explicitly expressed using Kirchhoff's formula, 

02 (t, x) = ——r / [tax (y) + a (y) + (y - x) ■ Va (y)] da (y) , 

4 ?rt JdB(x,t) 

where B(x, t) is the ball in M 3 of radius t centered at x, and where da is the standard 
surface measure for dB(x,t). In addition, from Duhamel's principle, 

47rJo JdB(x,r) r w JB{x,t) \V ~ x \ 

Note that this can be interpreted as an integral over the null cone segment 
Nq (t, x) = { (s, y) e R x K 3 | 0< s<t, \y-x\=t-s} . 

There exist many geometric analogues for the above model result, which provide 
possibilities for studying wave equations in other settings. An important early ex- 
ample is the result of S. Sobolev in [17] , in which a first-order parametrix was applied 
in order to prove well-posedness for second-order linear wave equations with vari- 
able coefficients. Y. Choquet-Bruhat made use of a similar representation formula 
in her celebrated local existence result for the Einstein- vacuum equations; see pQ. 
More recently, Chrusciel and Shatah, in [3] , used the Hadamard-type parametrix in 
[7] in order to extend the classical global existence result of Eardley and Moncrief 
for the Yang-Mills equations, cf. [H [5], to globally hyperbolic curved spacetimes. 

In contrast to the first-order variants of [U [17] , which we collectively refer to 
as "Kirchhoff- Sobolev- type" parametrices, the Hadamard-type parametrix of [7] 
achieved greater precision by making use of infinitely many derivatives of the metric 
g and requiring the convexity of the domain under consideration. For instance, 
since the geometric wave equation no longer obeys the strong Huygens principle, 
the Hadamard parametrix depends on all points in the causal past J~{z) of the 
point z under consideration. Consequently, geodesic convexity is required to make 
sense of the formula. This differs strongly from Minkowski spacetime, in which 
the fundamental solution in (1 + 3)-dimensions is supported entirely on the past 
null cone. These smoothness and convexity restrictions, however, can often be 
undesirable in the context of nonlinear problems in partial differential equations. 

We now focus our attention on [HI Thm. 3.11], which presented a Kirchhoff- 
Sobolcv parametrix directly handling covariant tensorial wave equations on arbi- 
trary (1 + 3)-dimensional Lorentzian manifolds. For convenience, we refer to this 
representation formula as KR. In particular, KR, applied at a point p G M, enjoys 
the following features: 

• KR directly handles the tensorial wave equation j2} in a completely co- 
variant fashion, in particular without reducing to a scalar wave equation. 

• The parametrix is supported entirely on the past null cone of p. Indeed, 
the formula is expressed as integrals of $ and ^ along this cone. 
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• Since the past null cone can degenerate away from p, then KR is valid only 
locally near p, on a "regular" portion VV~ (p) of the cone. In other words, 
the validity of KR is constrained by the regularity of the null exponential 
map, a weaker condition than the geodesic convexity assumption of [7|. 

• The formula KR contains "error terms" , expressed as integrals along Af~ (p) 
of $ along with various geometric quantities on J\f~ (p) . These are a result 
of the nontrivial geometry of the spacetime. 

• The formula KR depends only on quantities defined on J\f~(p) 1 i.e., it is 
independent of extensions of any quantities to neighborhoods of J\f~(p). 

• The formula KR can be systematically generalized to covariant wave equa- 
tions on arbitrary vector bundles over M with a bundle metric and a com- 
patible connection; see [12] Thm. 4.1]. 

For instance, [321 Thm. 4.1] can be used to reprove global existence for the 
Yang-Mills equations in R 3+1 in a similar manner as in [5]; see [12l Sec. 5]. In 
addition, since KR uses fully tensorial and covariant methods, then the above can 
be achieved without reference to Cronstrom gauges. This "gauge-invariant" method 
can likely be extended to treat curved spacetimes as well. 

Moreover, KR was applied in the proof of the breakdown criterion for the 
Einstein- vacuum equations in [14] , p The main result of |14j is the following: 

Theorem 1. Suppose {M,g) is an existing Einstein-vacuum spacetime, given as a 
constant mean curvature foliation 

M = |J S T) t <h< 0, 

i n <T<ii 

where each S r is a compact spacelike hypersurface of M satisfying the constant 
mean curvature condition tr k = r, and where k denotes the second fundamental 
form of ' E T in M. In addition, assume the following criterion: 

II fe II z,°°(M) + ll V ( lo g n )lli-(M) < °°- 
Then, (M,g) can be extended as a CMC foliation to some time t\ + e. 

In particular, KR was utilized to derive L°°-bounds for the curvature, which was 
shown to satisfy a covariant tensorial wave equation; see Q31 Sec. 5]. These L°°- 
bounds were essential for obtaining uniform energy estimates along the timeslices 
of the spacetime. Using local existence theory and various elliptic estimates, one 
could then derive Theorem [1] 

1.2. The Generalized Formula. The parametrix KR, however, is insufficient for 
handling the analogous breakdown problem for the Einstein-Maxwell (and similarly, 
the Einstein- Yang-Mills) setting. The main result, discussed in [16], is the following: 

Theorem 2. Suppose (M,g,F) is an existing Einstein- Maxwell spacetime, where 
F is a 2-form in M representing the Maxwell field. Assume M is given as a CMC 
foliation as in Theorem [3 and assume the following criterion: 

H fc llL~(M) + ll V (logn)|| L « (M) + II^IIl-(m) < °°- 
Then, (M, g, F) can be extended as a CMC foliation. 



The entirety of the proof of the breakdown/continuation argument stated in |14| . however, 
spans multiple papers, including [9l 1101 [TTl 1121 [T51 [Til [181 [19] . 
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The failure of KR in the proof of Theorem [2] is due to the presence of nontrivial 
first-order terms in the wave equations for the curvature and the Maxwell tensors. 
More specifically, if we apply KR, then we obtain null cone integrals containing 
"bad" components of the derivative of the curvature, which we cannot control a 
priori using standard local energy estimates. We also obtain equally untreatable 
"bad" components for the Maxwell tensor. 

In order to obtain Theorem [2] we must deal with the aforementioned first-order 
terms. To achieve this, we generalize KR to handle a system of tensor wave equa- 
tions of the form ((T|). In particular, we must handle the first-order portion of (TT]) 
differently. As we shall see later, the main "trick" will be to eliminate those terms 
containing derivatives of the $ (m) 's which are transverse to the null cone. This is 
intuitively unsurprising, since such derivatives, in both flat and curved spacetimes, 
are generally the most difficult to control using local energy estimates. 

Before delving fully into technical details, we first compare the approach of this 
generalized representation formula to that of its predecessor KR: 

• Similar to developments in [7], we handle the contributions of the first- 
order terms by suitably altering the transport equation in [T^j. We add 
corresponding terms to these transport equations, which include the P {mc ) 's. 

• Due to coupling, all the equations in (jTJ must be handled concurrently. 
Consequently, the transport equation of [12] becomes a coupled system of 
n transport equations. 

• Like for KR, the parametrix of this paper will be supported on regular past 
null cones and will depend entirely on quantities defined on these cones. 

We opt for a different approach to the derivation of our formula than for that of 
KR. The main differences in principle between the two proofs are as follows: 

• A major component in deriving KR is the optical function, whose level sets 
are null cones. Here, we avoid any reference to such a function. As a result, 
we can eliminate some extraneous geometric assumptions needed in |12j . 

• While pj] makes heavy use of distributions on manifolds in an informal 
fashion, we take a more technically rigorous path and remain at the level 
of the calculus operations represented by such distributions. In particular, 
rather than dealing with derivatives of the ^-distribution like in |12) . we 
deal with corresponding integrations by parts. 

• In our derivation, we integrate by parts only the derivatives in directions 
tangent to the null cone, while in [12], integration by parts is applied lib- 
erally to derivatives in all directions. The latter approach results in terms 
depending on quantities off the null cone, which must then be cancelled 
through additional meticulous integrations by parts. By exercising further 
discretion, we can avoid such troubles altogether. 

• In this paper, we will make explicit the dependence of the parametrix on 
the initial data. □ This is in contrast to [12], which focused only on the 
case of vanishing initial data. 

As mentioned before, we obtain subtle improvements over the assumptions re- 
quired for KR to be valid. Assumptions Al and A2 of [12] Sec. 2.1], assumed 
when deriving KR at a point p € M, postulated that certain local hyperbolicity 



By "initial data", we mean the values of the <J>< m )'s and D<I>< m >'s on a spherical cross-section 
of the past null cone. These correspond ag and a\ in the scalar Minkowski model problem. 
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and null regularity conditions hold for all points in some neighborhood of p. In 
contrast, in our derivation, since we avoid the optical function and work exclusively 
on the null cone, we require only the null regularity condition for p. 

We also previously noted that both the representation formula of this paper and 
KR depend only on quantities defined on the null cone. Since the derivation of 
KR introduces a multitude of terms which depend on quantities off the cone, that 
the final result does not depend on quantities off the cone appears as a "miracle" 
resulting from numerous cancellations. In contrast, in our derivation, we automati- 
cally have at each step of the proof that every term is dependent only on quantities 
defined on the null cone. Therefore, our derivation highlights this property of the 
parametrix as a natural rather than surprising consequence. 

1.3. The Main Result. Now, we state an abridged schematic version of the main 
theorem of this paper. The precise statement will be deferred until Theorem[7J after 
we develop the background and notations needed to fully describe the parametrix. 

Theorem 3. Consider the system (flj on (M,g), let W denote a "regular" segment 
of the past null cone about p G M , and let J (1) , . . . , J (n) be tensors at p of the same 
ranks as $ (1) , . . . , $ ( "' ; respectively. Then, we have the schematic formula 

n n p 

y^ J (m) • $ (m) | p = -^2 ^ <m> ' * <m) + ( error terms) + (initial data) , 

where: 

• The A (m) 's are tensor fields on J\f which satisfy a system of transport equa- 
tions depending on the J (m) 's and on the geometry of N ' . 

• The error terms can be expressed as integrals of quantities on M ■ 

• The initial data contribution can be expressed as integrals over the lower 
boundary of Af (i.e., a spherical cross-section of the past null cone). 

In addition, in the final section, we will discuss extensions of the parametrix to 
covariant wave equations on sections of vector bundles, similar to |12[ Thm. 4.1], 
but once again with nontrivial first-order terms. We also discuss the main theorem 
in the context of this extended framework, and we discuss how this extension can 
be applied to treat the Einstein- Yang-Mills breakdown problem. 

Acknowledgements. The author would like to thank Professor Sergiu Klainer- 
man for his insights regarding this topic, and for hours of helpful discussions. 

2. Geometry of Regular Null Cones 

Fix p £ M, and consider the null exponential map exp about p, defined as the 
restriction of the exponential map about p to the past null cone 91 of the tangent 
space T p M. [] Define the past null cone N~ (p) of p to be the image of exp . Note 

that N~(p) is ruled by the past inextendible null geodesies beginning at p. 

We call a neighborhood U in 91 of the origin star-shaped iff for each q € U, the 
open segment between the origin and q is also contained in U. Fix such a star- 
shaped U C 91, and assume exp is also a global diffcomorphism on U . From now 



Here, we assume *Tt does not include the origin of T p M. 
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on, we shall denote the image exp (U) by Af~(p), which we refer to as the regular 

past null cone, u In particular, M~ (p) is a smooth null hypersurface of M. 

Remark. Recall that since one can find convex subsets about any point, then one 
can always construct such a set J\f~ (p) . 

Remark. For nonpathological spacetimes, we can systematically construct such 
Af~(p) via global considerations. For instance, in the globally hyperbolic setting, 
we can determine A/" - (p) from the past null boundary J~ (p) — I~ (p) , where I~ (p) 
and J~ (p) are the chronological and causal pasts of p, respectively; see [H 115) . In 
practice, this is the situation adopted in most applications. 

We remark here that a loss of such regularity occurs at what are called terminal 
points of N~(p). At such a point z, one of the following scenarios hold: 

• Two distinct past null geodesies from p intersect at z. In other words, the 
map exp fails to be one-to-one at z. 

• The pair p and z are null conjugate points. In other words, the map exp 
fails to be nonsingular at z. 

From now on, we shall deal only with the regular portion A/" _ (p), on which 
we have a smooth structure. Furthermore, since tangent null vectors in J\f~(p) 
have vanishing Lorentzian "length" and are orthogonal to Af~(p), they cannot be 
normalized without introducing vectors transversal to J\f~(p). Consequently, in 
our treatment, we will require an additional choice of a future timelike unit vector 
t G T P M. Here, we fix such a t for the remainder of this paper. 

2.1. Null Generators. We define the null generators of N~(p) to be the inex- 
tendible past null geodesies 7 on M which satisfy 7(0) = p and gr(7'(0),t) = 1. 
We can smoothly parametrize these generators by S 2 using the following process. 
If we choose an orthonormal basis eo, . . . , e^ of T p M, with eo = t, then we can 
identify each w G S 2 with the null generator 7^ satisfying 7^,(0) = — e +uj k ek- For 
convenience, we assume such a parametrization of the null generators of A/" - (p) , 
and we denote by j u , ui G § 2 , the null generator corresponding to u. 

Remark. The objects on M~{p) that we will discuss can of course be defined in- 
dependently of any parametrization of the null generators. However, for ease of 
notation, we will work explicitly with § 2 . 

In addition, we define the vector field L on A/"~ (p) to be the tangent vector fields 
of the null generators, i.e., we define L\ loj ^ — 7^,(i>) for any u G § 2 . Note in 
particular that L is a geodesic vector field. 

2.2. Spherical Foliating Functions. We say that a function / G C ca (N~{p)) is 

a foliating function of Af~(p) iff / satisfies the following conditions: 

• For every to G § 2 , 

li x m n /k(.)=°- 

• The function 

(4) $ = $ f = (Lfy 1 G C°° (AT- (p)) 

is everywhere positive on Af~(p). In other words, / is strictly increasing 
along every null generator of J\f~ (p) . 



"In other words, we can think of M (p) as a "past null convex" set. 
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limtfL ^ =^o- 



There is a constant do > such that for all uj G § 2 

• For every w£§ 2 , the limit 

exists and is finite. 
The function ■& defined in ((4]) is called the null lapse of /. 

Remark. We can in fact weaken the third condition in the above definition so that 
•& converges to different limits along various null generators. This generalization, 
however, is currently without practical applications and complicates various initial 
limit computations. Therefore, we restrict ourselves to the simpler definition. 

In the remainder of this paper, we will let / denote an arbitrary "default" fo- 
liating function for Af~(p), and we let § = $/. The positivity of •& implies df is 
everywhere nonvanishing, so the level sets of /, denoted 

S v = Si = {qeM- (p) \f(q)=v}, v>0, 

form a family of hypersurfaces of Af~(p). Since L represents the unique null direc- 
tion tangent to M~ (p) , and the positivity of $ implies that L is transverse to each 
S v , we can conclude that each S v is spacelike, i.e., Riemannian. We shall adopt the 
following notational conventions: we let A and V denote the induced metrics and 
the Levi-Civita connections on the <S„'s, respectively. 

Define the quantity i(p) = i/(p) to be the supremum of all v > such that for 
every u e § 2 , there exists a point q on both 7^ and 7V _ (p) such that f{q) = v. 
This can be interpreted as N~(p) remaining regular up to /-value i(p). \J Note that 
the definition of i(p) and the positivity of t? imply that S v is diffcomorphic to S 2 
whenever < v < i(p). 

In addition, for any < Vi < i(p), we define the null cone segment 

Af- (p; t>i) = Nj (p; «i) = {q G N~ (p) \ f (q) < Vi} , 

while for any < v\ < V2 < i(p), we define 

M~ (p; vi,v 2 ) = Mj (p; vi,v 2 ) = {g£^ (p) | v x </(<?)< v 2 } . 

The most natural example of a foliating function is the affine parameter s on 
7V~(p), given by s(7 w (w)) = v for each w G § 2 . Note in particular that s converges 
to at p and i? s = 1 everywhere. The foliation of Af~ (p) by the <S^ 's is called the 
geodesic foliation. 

Remark. Note that s depends on the normalization t. 

If the spacetime M is foliated by a global time function t, then we have another 
natural foliating function t p for 7V _ (p), given by t p (q) — t(p) — t{q). This is the 
foliating function used in the breakdown results of [331 US] . 



If M (p) is defined in terms of the global geometry of M, as was discussed before, then i(p) 
corresponds to the "null injectivity radius" of M~ (p) with respect to /. 
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2.3. Tensors. Here, we list the notations we will use to describe various tensor 
fields on M~{p). Again, we assume an arbitrary foliating function / for J\f~(p). 
We begin with the following objects: 

• A tensor w at some q € S v is said to be horizontal iff w is tangent to S v . 

• We denote by T ~J\f~ (p) the horizontal bundle over Af~ (p) of all horizontal 
tensors of total rank k at every q 6 J\f~(p). 

• Similarly, we denote by T Af~ (p) the extrinsic bundle over Af~ (p) of all 
tensors in M of total rank I at every q £ ]\[~{p). 

For a vector bundle V, we let TV denote the space of all smooth sections of 
V. By this formalism, then TT k J\f~(p) and TT N~{p) denote the spaces of all 
horizontal tensor fields of rank k and extrinsic tensor fields of rank I on M~ (p) , 
respectively. In other words, a horizontal tensor field A G TT_ Af~(p) smoothly 
maps each q 6 J\f~ (p) to a horizontal tensor of rank k at q; extrinsic tensor fields 
can be characterized analogously. For example, the restrictions to N~(p) of the 

2 

spacetime metric g and the curvature R can be treated as elements of TT Af~ (p) 
and TT Af~(p), while the induced metrics A on the iS„'s can be treated as an 
element of TT 2 J\f-{p). Lastly, note that TT°Af-(p) = TT°Af-(p) = C°°(Af-{p)). 
Next, we define the mixed bundles over M~ (p) to be the tensor product bundles 

T k T l N- (p) = T k N- (p) ® T l Af- 0) . 

Similarly, we will call a section A e TT_ T J\f~ (p) a mixed tensor field on A/"~ (p) . 
By the duality formulation, we can consider such a field A as a bilinear map 

A : rT fe AT (p) x IT W~ (p) -» C°° (M~ (p)) , 

or as the corresponding C°° {N~ {p))~ valued multilinear map on k horizontal vector 
fields and I extrinsic vector fields. 

In index notation, we adopt the following conventions: 

• Horizontal indices will be denoted using Latin letters and will take values 
between 1 and 2, inclusive. 

• Extrinsic indices will be denoted using Greek letters and will take values 
between 1 and 4, inclusive. 

• Collections of extrinsic indices will be denoted using capital Latin letters. 

2.4. Covariant Differentiation. Recall that the Levi-Civita connection D on M 
induces a connection D on the extrinsic bundles T Af~ (p) . For any A e TT Af~ (p) 
and a vector field X on Af~(p), we can arbitrarily extend A to a neighborhood of 
A/"~ (p) and define DxA = DxA. It is easy to see that this definition is independent 
of the chosen extension of A. Moreover, it is clear from definition that Dxg = 
for any vector field X on A/"~ (p) . u 

The Levi-Civita connections y on the <S,/s naturally aggregate to define a con- 
nection V on the horizontal bundles TT~J\f~(p). Let Xka vector field on Af~(p): 

• For a scalars / € C°°(J\f~(p)), we define y*x/ = Xf, as usual. 

• For a horizontal vector field Y, we define y X Y to be the projection onto 
the <S„'s of Z?x^- Note that if X is horizontal, then y '^7 is precisely the 
corresponding covariant derivative with respect to the iS„'s. 



"Technically, by Dxg, we mean Dx acting on the restriction of g to J\f (p). 
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• For a fully covariant A E FT k Af-(p), then f x A E TT k J\f-(p) is naturally 
given as follows: for horizontal vector fields Y±, . . . , Y k , 

f x A (Fi, ...,Y k )=X[A(Y l ,..., Y k )\ - A (f x Y u Y 2 , ..., Y k ) 

-...-A(Y 1 ,...,f x Y k ). 

This definition of V generalizes the usual ^-covariant derivative on the S v 's to also 
include the L-direction. Note also that V X X = for any vector field X on Af~(p), 
i.e., V remains compatible with the horizontal metrics. 

We can canonically combine the connections y and D on the horizontal and 
extrinsic bundles to obtain generalized connections y on the mixed bundles. The 
basic idea is to have V behave "like y on the horizontal components" and "like 
D on the extrinsic components" . More specifically, if A E TT_ T A/ -- (p) is a fully 
covariant mixed tensor field, X is a vector field on Af~(p), Yi,. ,., Y k E FT}N~(p) 
are horizontal vector fields, and Z\,...,Zi E TT M~{jp) are extrinsic vector fields, 
then we define ^/ x A £ YT k T A/ -- (p), interpreted as a multilinear map, by assigning 
to the expression *$ X A(Y\, , . . ,Y k ; Z\, . . . , Z{) the following value: 

X [A (Y u ...,Y k ;Z 1 ,.. .,Zi)\ - A (H X Y X ,Y 2 , ...,Y k -Z u ...Z t ) - ... 

-A(Y X ,..., Y k _ x ,y x Y k ; Z 1 ,...,Z l )-A(Y 1 ,.. .,Y k ;D x Z 1 ,Z 2 , ...Zi) 

-...-A(Y 1 ,...,Y k ;Z 1 ,... 7 Z^ u D x Zi) . 
The main observations concerning the above construction are the following: 

• The mixed covariant derivatives satisfy Leibniz rules similar to ^ and D- 
covariant derivatives. 

• For any vector field X on 7V _ (p), both V^ A = and V x9 — 0- 

The primary consequence of these observations is that the same integrations by 
parts operations can be justified for ^-derivatives of mixed tensor fields as for V- 
derivatives on horizontal tensor fields. This property was implicitly used in |121I14) . 
Finally, we will make use of the following notations: 

• For any A E TT k Af~{p), we let fA E TT k+1 Af-(p) denote the horizontal 
tensor field mapping a horizontal vector field X to ^ X A. This is precisely 
the covariant differentials of A on the S^'s. 

• Similarly, for any A E TT k T l J\f-{p), we let fA E TT k+1 T l JV_{p) denote 
the mixed tensor field mapping a horizontal vector field X to ^fi X A. 

• We also define the horizontal and mixed Laplacians in the usual way: 

& = \ ab f ab , J=X ab f ab . 

• Consider also the operators y j = d^ ' L and ^ f = -O^l- These are horizon- 
tal and mixed covariant derivatives in the tangent null direction, subject to 
the normalizations ^// = y // = 1- 

For further details involving the above constructions, see [T6j Sec. 1.2]. 

2.5. Parametrizations. We can parametrize J\f~(p) by the foliating function / 
and a value ui E S 2 . For < v < i(p) and lu E S 2 , we can identify (v,uj) with 
the unique point q on the corresponding null generator j u with f(q) = v. As a 
result, we can naturally treat any (f> E C co (M~(p)) as a function of / and w. For 
any such (j>, we denote by 4>\(v,u>) the value of </> at the point q corresponding to the 
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parameters (v, lu). We will freely use this (v, w)-notation throughout future sections 
without further elaboration. 

2.6. Null Frames. In general, null frames are local frames I, m, ei, e^ which satisfy 

9 [t, I) = g (m, m) = 0, # (l, mj = -2, 

5 M~> e aJ = 3 (»w 5 e a) = 0, 3 (e a , e fc ) = <5 ab , 

Here, we define null frames which are adapted to the /-foliation of Af~(p). 

Each point of S v is normal to exactly two null directions, one of which is repre- 
sented by L. We define L, called the conjugate null vector field, to be the vector 
field in the other normal null direction, subject to the normalization g(L,L) = —2. 
Next, we append to L and L a local orthonormal frame ei,e2 on the <S„'s. Then, 
{L,L, ei, ei\ defines a natural null frame for M~ (p). 

In this paper, we will index only with respect to adapted null frames. 

• Horizontal indices 1, 2 correspond to the directions e\ and ei- 

• L corresponds to the index 3, while L corresponds to the index 4. 

2.7. Ricci Coefficients. We will make use of the following connection quantities: 

• Define the null second fundamental forms x, X £ rT 2 7V~ (p) by 

X (X, Y) = g (D X L, Y) , x (X, Y) = g (~D X L, Y) 

for any horizontal vector fields X and Y on Af~ (p) . Note that x an( i X are 
symmetric, since both L and L are normal to the S v 's. 

• We often decompose x m to its trace and traceless parts: 

tT X = X ab Xab, X = X-7;( tT x) X - 
We will also use an analogous decomposition for x- 

• Define C,V € TT^N'fr) by 

C (X) = lg (D X L, L) , r L (X)= ig (X, D L L) 

for any horizontal vector field X on Af~(p). 
The quantities tr^, X\ and £ are called the expansion, shear, and torsion of Af~(p) 
(with respect to the /-foliation), respectively. 

In addition, we have the following relation between £ and rj: 

(5) r [ =-( + f(logtf). 
For a proof, see [9j Prop. 2.7]. 

Given a null frame L,L, e±, e2, a vector field Z G TT Af~(p) is decomposed 

1 1 2 

£ = -^{ZiDL- -g(Z,L)L + ^2g(Z,e a )e a . 

Using this identity, we can decompose covariant derivatives along Af~(p): 

(6) D a e b = f a e b + -^X + -^XabL D a L = Xa&b - QL, 

D a L = x a fc e b + Co£, D L e a = f L e a + i I( L, 
D L L = 0, D L L = 2r ] a e a . 
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Lastly, we define the mass aspect function fj, £ C°° (J\f~ (p)) by 

(7) H = U\a \T b X ab + ICI 2 + "#4343 - ^43, 

and we note the following transport equation satisfied by tr \: 

(8) t L (tr%) = 2f\ + 2 ^l 2 - i (tr X ) (tr%) - X ab X ab + ^4343 - #43- 
For details on ©, see g] or QU Sec. 3.2]. 

2.8. Integration. Since Af~(p) is null, we have no volume form on Af~(p) with 
respect to which we can integrate scalar functions. However, we can still give a 
natural definition for integrals of functions over J\f~(p), as long as we have a fixed 
normalization for M~{p). Indeed, we define this integral by 




(9) 

IN-{ P ) 

for any </> £ C°°(N~(p)) for which the right hand side is well-defined, where the 
S^ 's are the level sets of s on M~ (p) . We can similarly define integrals over any 
open subset of Af~(p), in particular the segments M~{p]V2) and Af~(p;Vx, v<z). 
By a change of variables, we can restate © in terms of any foliating function /: 

Proposition 4. For any foliating function f of J\f~(p), we have 

tf-Adv 



/AA-(p) J a \Js v J 

for any integrable 4> £ C°° (J\f~ (p)) . 

We will also need the following derivative formula. Here, the integrals over S v 
are the usual integrals over Riemannian manifolds. 

Proposition 5. Let <f> £ C° c {N^{p)). Then, for < v < \{p), 
j- f <f> = f [ff<t> + #(tr X )<f>]. 

For proofs of the above propositions, see [16l Sec. 3.4] and [9]. 

Remark. We can justify the definition @ using energy estimate considerations. 
Suppose V is a region in M , with a portion of its boundary given by a neighborhood 
V' in Af~(p). Then, if X is a 1-form on M, and we integrate D a X a over V and 
apply the divergence theorem, then the boundary integral corresponding to V' is 
precisely given by 

(10) I g(X,L). 

Jy 

This formula is fundamental to the local energy estimates applied in [141 116] . 

Remark. Note that the definition @ depends on the normalization t chosen for 
Af~(p), while the argument in the previous remark establishes that the expression 
(|10|) is in fact independent of t. 
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2.9. Initial Values. Another important set of properties we shall need in our 
analysis concerns the "initial values" of many of the objects we have defined above, 
that is, we wish to compute their limits at p along the null generators of N~(p). 
We list the general results for arbitrary foliations below. Their proofs involve 
applications of convex geometry among other technicalities, hence we omit them in 
this paper. For details, consult [161 Sec. 3.3]. An earlier account for the special 
case of the geodesic foliation is presented in p~8] . 

Proposition 6. The following limits hold for each u> £ E> 2 . 

• We have the following comparisons between f and the afjine parameter s: 

Here, i?o is the initial value of the null lapse #. 

• For each integer k > 0, 



(11) 



lim - 

v\a f 



(v,ui) 



(12) 



lim v k 1 

v\0 



(v,ui) 



= 0, 



lim v k 



f k $ 



(v,u>) 



= 0, 



(13) 
(14) 
(15) 



where ^/ denotes the k-th order horizontal covariant differential, i.e. 
operator rf applied successively k times, [j 

• We have the following limits for \- 

lim Ii9(trx) -2/- 1 ||, . = 0, lim |x|L, ^ = 0. 

• We have the following limits for £ and rj: 

l™-ICII ( ^-limH!Zll^)=0- 

• We also have the following limits for \ : 



the 



limw(trx)L , =-2tfr 1 J 



lim v v , . = . 



v\0 v — ' ' ( 

Let <f> G C°°(Af-(p)), let O G C°°(§ 2 ), and suppose 

w'e§ 2 . 



lim 

v\0 



\(v,u') - <W\(u')> 
Then, the following integral limit holds: 



?9 2 



Mu,'^ 1 ' 



(16) lim v~ 2 

"V js v Js 2 

where the integrals on the right hand sides are with respect to the standard 
Euclidean measure on S 2 . 

In the previous proposition, the tensor norms | • | denote the natural Ricmannian 
norms on the 5,,'s, that is, for any A £ TT_ Wf~(p), we define 

\A\ 2 =X(A,A) = (f[X a A A ai ... ak A bl ... bk . 

Remark. We can also see from Proposition\^that the behaviors of the objects listed 
within the proposition all "tend to their corresponding values in Minkowski space" 
at p. For example, the expansion tr \ in Minkowski space is precisely 2s _1 , which 
is the asymptotic behavior near p demonstrated in (Q2 



"Note that V 7 here is defined with respect to the /-foliation. 
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Remark. We can also derive initial limits for derivatives of the Ricci coefficients. 
For details, see |16[ Sec. 3.3]. 

Remark. The initial values for the Ricci coefficients behave even better than in 
Proposition^ in the case of the geodesic foliation; see |16[ Sec. 3.3] or |18j . 

3. The Main Theorem Statement 

Before we can state the representation theorem, we must make a few additional 
definitions. Again, we suppose Af~ (p) is normalized by t £ T p M and foliated by 
/ G C ca (N^(p)). Recall that from t and /, we can define the null vector fields L 
and L, as well as local null frames L, L, ei, e-i adapted to the /-foliation of Af~(p). 

For each m £ {1, . . . , n}, we define the extrinsic tensor fields 

B im) eTT r{m> M-{p) 

such that they satisfy the coupled system of transport equations 

(17) f f B (m)I = -- 0(tr X )-i B (m)J + -^P (om) 4j 7 B (o)J , l<m<n 



f 



2 

c=l 



along the null generators of A/" (p), where I and J here denote collections of r (m) 
and r (c) extrinsic indices. We also stipulate the initial conditions 

(18) B im) \ p = J {m \ l<m<n, 

where each J (m) is an aribtrary tensor of rank r (m) at p. Note that the validity of 
the system (jTTJ), (fTS]) follows from the initial limit (|TU| for tr%. 
We also define the extrinsic fields 

a™ = f-^B™ g rr r(m> Ar ( P ) , 

and we note that the A {m) 's satisfy the transport equations 

— 1 1 - 

(19) f L A^ I = --(tr X )A^ I + -Y / P (cm> ij I A^ J , l<m<n. 

c=l 

The A (m) 's are the analogues of the tensor field A in [12] and correspond to 
the factor r _1 = \y — x\~ l in the explicit expression for <j>i in the model problem. 
Moreover, if n = 1 and the P< mc >'s vanish, then (fT9|) reduces to the transport 
equation given in [12]. Note in particular that even if a specific J (m) vanishes, v4 (m) 
can still be nontrivial due to the coupling in (|19l) . 

Lastly, we define the coefficients 

z/ cm) e TI^ C +r m Af~ (p) , l<m,c<n 
by the formula 

(20) i/"»y = -Tp' cm \/ + ^y 4 -p Com, 3/ + cp^a/ + \ (trx) p'^^j 7 



+ \ (tr X ) P (cm) 3 / + ^ E /J(rf, 4/ rJ,i "' , 3/- 

d=l 

These tensor fields will be present in the error terms of the representation formula. 
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3.1. The Main Theorem. We are now ready to state the main theorem. 

Theorem 7. Assume the following: 

• Let n be a positive integer, let r (1) , . . . , r (7l) be nonnegative integers, and 
suppose for each 1 < m, c < n, we have defined tensor fields <j> (m) ; \p (m ' ; 
and P (mc) on M of ranks r (m> , r <m) , and 1 + r (m) + r (c) , respectively. 

• Suppose the <j> (m) 's, \I/ (m ' 's, and P< mc > 's satisfy the system JlJ. 

• Fix p £ M, and suppose a regular portion A/" - (p) o/ f/ie past null cone 
N~(p) is normalized and foliated by t £ T p M and f 6 C°° (Af~ (p)) . 

• Let vq be a constant satisfying < Vq < i(p). 

• For each 1 < m < n, we define the extrinsic tensor fields 

£(m> e rT r<m V^ (p) , A'- m) = f-^-B™ e rr r<m V" (p) , 

along with a tensor J (m) ofrankr (m) atp, such that the systems of transport 
equations dTTJ), (H]), dHJ) ZioZd. 

Then, we have the representation formula 

(21) 47n? ^ j<'" )7 $^j| p = 5 (p; «o) + e 1 (p; «o) + e 2 (p; «o) + 3 (p; wo) , 

m— 1 

tf/iere 

• T7ie "fundamental solution term" $(p;vo) is given by 



(22) 






• 



TTie "principal error terms" € (p;i>o) Are given by 



IT, r. 

(23) e 1 (p; vo) = -J2 fA'-'^Wi 

m=1 W-(p;t>o) 
n « 

+ E/ (c°-n a )^«>4 Cm)/ * <m) i- 

• TTie remaining "error terms" <£ 2 (p;vo) are given by 

(24) £ 2 (p;i; )=^ / M ■ A'"" 7 ^, 

n . 

rn,c=l"' JV ' _ (P;' u o) 

n r. 

m,c=l^ _ (Pi'"o) 

where i?43[$ (m) ]/ denotes the curvature quantity 

i? 43 [* <m> ] 7 - ^43* <m, 7 - ^34$ (m, / ; 

and where "D43 " and "D34 " denote second covariant derivatives. 
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• The "initial value terms" 3(p;vo) are given by 

(25) 3 (p; v ) = ~ J2 J ( tr ^) A (m)J <I>« 7 - J2 I ^ ! D^ 

i " r 



(m, 7 



C=1 J5„ 

Here, we have indexed with respect to arbitrary null frames L.L, ei,e2 adapted to 
the f -foliation. The capital letters I, J refer to collections of extrinsic indices. 
Furthermore, the error terms (£ (p;i>o) can oe alternately expressed as 

n „ n r. 

(26) £ 1 (p;v )=Y, JA^ I ^^ I + 2J2 C -^ a A (m) T <$> {m) i . 

Note that §(p; vo) corresponds to the explicit form of <pi in the model problem, 
while (£ (p;i>o) and (£ 2 (p;uo) are error terms which vanish in the model problem. 
Lastly, 3{p;vo) describes the initial value contributions, which correspond to the 
KirchhofF formula for <f>2 in the model problem. 

Remark. Although the representation formula was stated in (|21 [) - f|26[) in index 
notation, this was done only as a matter of convenience. It is easy to see that these 
expressions can in fact be described invariantly. 

Remark. In particular, we can use (|21[) to examine the value of any <j> (m) \ p indi- 
vidually by setting J <c) = for all c =/= m. 

3.2. Improvements over KR. Theorem [7] contains a number of improvements 
over the parametrix KR. The most apparent is that Theorem [7] handles the system 
dl]) and not just the single equation ([2]). While [12] only presented the case with 
vanishing initial data, here we present the general case for arbitrary initial data; 
the dependence on the initial data is given explicitly by the J(p; ^o)-term (|25j) . 

Remark. In fact, the extended parametrix presented in [121 Thm. 4.1] can handle 
the "multiple wave equations" aspect of {T]) 7 but it does not cover the issues of 
general first- order terms. For details, see section 5. 

Additionally, Theorem[7]weakens the assumptions required for the representation 
formula to be valid. The main assumptions needed for KR are the local spacetime 
conditions Al and A2 given in [121 Sec. 2.1]. First, since we make no reference to 
the optical function and remain entirely on A/" - (p) in our derivation, we only need 
that the null cone regularity stated in A2 hold for A/" - (p) , rather than for A/"~ (q) 
for every q in a neighborhood of p, as needed in |12) . This assumption is reflected 
in Theorem [7] by working only within the segment M~(p;vo), where i>o < i(p). 
Moreover, since we work exclusively on A/" - (p), the local hyperbolicity condition 
Al is entirely superfluous here. 

3.3. Exclusive Dependence on the Null Cone. We also previously noted that 
the representation formula depends only on quantities defined on J\f~(p). Indeed, 
by inspecting the individual portions (|22 |) - (|24|) . we can immediately see that all the 
objects within the integrands can be expressed as horizontal, extrinsic, or mixed 
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tensor fields on A/" - (p) . The representation formula remains unaffected by potential 
extensions of these quantities off J\f~(p). lj 

We remark that this property of the parametrix does not conflict with the fact 
that wave equations on curved spacetimes do not satisfy the strong Huygens prin- 
ciple. This is demonstrated in Theorem [Jj by the recursive error terms (|23[) , (1241) , 
(f26l) . which contain the unknowns $ (1) , . . . , $ <7l) themselves. 



3.4. Extensions Beyond Cut Locus Terminal Points. As stated in Theorem 
[7J the parametrix is valid only when the past null exponential exp remains a 
global diffeomorphism. We can, however, trivially extend Theorem [7J beyond cut 
locus terminal points, that is, terminal points resulting from the intersection of two 
distinct null generators, by lifting the representation formula to the tangent space 
T p M via exp . By pulling back over exp , we can consider all the tensorial objects 
referenced in Theorem [7J as fields on the past null cone 9T of T p M. At the level 
of the tangent space, cut locus terminal points no longer exist, so we can apply 
Theorem [Jj as usual to obtain a representation formula on 91. 

Using exp , we can express the above representation formula on 9t back in terms 
of the past null cone of p, beyond any cut locus terminal points. Due to intersecting 
null geodesies, some objects in the integrands, such as the Ricci coefficients, will 
be multi-valued at cut locus points. However, in this extended formulation of 
Theorem [7J we longer have the standard local energy estimates which were crucial 
to the general relativity applications [T3J [TB] . 

On the other hand, conjugate terminal points on past null cones pose a much 
more serious problem. Such points by definition are preserved by the lifting via 
exp to the tangent space T p M . Moreover, conjugate points are accompanied by 

the degeneration of the expansion tr \ and hence the A im) 's. Consequently, we have 
no natural extension of Theorem [7] which survives beyond conjugate points. 

3.5. Applications to General Relativity. Here, we briefly describe how Theo- 
rem[7Jis applied to the breakdown problem for the Einstein-Maxwell equations, i.e., 
to Theorem [5] First, in the vacuum analogue, one crucial estimate is the uniform 
bound for the spacetimc Ricmann curvature R, which satisfies a wave equation 

(27) a g R^R-R. 

The parametrix KR was applied to (J2JJ) in order to control R at a point p; the 
principal term to be estimated as a result is an integral over J\f~ (p) of the quadratic 
nonlinearity R ■ R. £j This was handled using a similar strategy as in the classical 
work [5] of Eardley and Moncrief, in that at least one of the R's in the integrand 
could be controlled by the flux density of R on Af" (p) . 

The Einstein-Maxwell case is similar, except that we now have a system of two 
wave equations for both R and the Maxwell 2-form F on M: 

(28) U g R = F ■ D 2 F + (R + DF) ■ (R + DF) + l.o., 

D g DF = F-DR+(R + DF) ■ (R + DF) + l.o.. 

We can apply Theorem [JJ to this system, and we can handle the quadratic non- 
linearities in the same manner as in the vacuum analogue. The first-order terms 



The only exceptions to these rules are the factors D^ 1 -" 1 ^ in the "initial data" terms 1 )251 . 
which depend on transversal derivatives of the <E>' m ''s. 

The analogue A of the transport equation in KR also appears here. 
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F ■ D 2 F and F ■ DR present a new challenge, however, as the methods in [M] us- 
ing KR fail here. However, these terms can be handled directly using Theorem [3 
which absorbs the effect of these terms into the transport equation. 

More specifically, in the setting of Theorem[71 we take $ (1) = R and $ <2) = DF. 
From (|28|) . we see that the first-order coefficients P (11) and P {22} vanish, while the 
coefficients P {12) and P i21) are schematically of the form F. Due to the extra control 
present for F stipulated in Theorem [21 the A {rnV s, which depend on the P (cm) 's, 
can still be adequately controlled. For details, see [IB] . 

4. Proof of the Main Theorem 

In this section, we assume the hypotheses of Theorem [3 and we derive in detail 
the parametrix of Theorem [7] Let < e < vq, and define for convenience 

K = Af- (p; e, v ) = {q e Af~ (p) \ e < f (q) < v Q } . 

4.1. General Integration Formulas. The most fundamental steps in the deriva- 
tion of Theorem [7] involve integrations by parts. Here, we state the general lemmas 
representing these steps. These lemmas deal with two separate cases: integrations 
by parts involving y , and integrations by parts involving tf L . 

We first examine the horizontal case. This is a simple consequence of the com- 
patibility of the mixed covariant differential V with both the spacetime metric g 
and the horizontal metrics A. 

Lemma 8. For any integer r > 0, S £ TT J\f~(p), and T G TT_ T N~{p), 

f S I -TT aI = -f Ts'-Ta!- [ ^(log^-S'T^, 

JM C JN e JM t 

where I denotes a collection of r arbitrary extrinsic indices. 
Proof. Applying Proposition U yields 

(29) j S 1 ■ f a T aI = J° (J # ■ sTtJ\ dv 

f T^-S'Taj) dv- j" (j #V -T a !\dv 

~ r '(/ ^ a ^^ ' siTai ) dv ' 

In the last step, we used that both V 3 and ^7 vanish. By Proposition [4] again, 
the last two terms on the right-hand side of (|29p are precisely the desired terms. 
Defining the horizontal 1-form to by uj(X) — § ■ S I T a jX a , then by the divergence 
theorem, the first term on the right hand side of (|2"9"|) becomes 

f a uj a ] dv = 0. D 

As a special case of Lemma 13 we have the following: 
Corollary 9. For any integer r > and S,T E TT Af~(p), 

f S 1 ■ jTi = - f TS 1 ■ f a Ti - f f a (log 1?) • S'faTj, 
JN e JM, JM C 

where I denotes a collection of r arbitrary extrinsic indices. 
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Next, we give the result for the null direction case: 
Lemma 10. For any integer r > and S, T 6 TT Af~(p), 



S 1 ■ f L T t 



UlS 1 ■ Tj 



(tr X )-S 7 T/- / S T Ti 



S'Tj, 



where I denotes a collection of r arbitrary extrinsic indices. 
Proof, lj First, we have 

/ S 1 -f L Tj = [ f L (S'Tj) - f t L S ! • TV. 

JAf e JM, JAf e 

By Propositions 2] and [5J the first term on the right-hand side can be handled by 



a/; 



Ul (S'Ti) 



"o 



S v 



f f (S'Tj) 



dv 



°^{! S S-T, 



dv' 



dv 



i)(tr X )S I T I 



s v 



dv, 



where we recalled that V / denotes the covariant derivative $ t l . The first term on 
the right-hand side is precisely the desired "boundary terms" . Finally, an applica- 
tion of Proposition @] to the second term completes the proof. □ 

4.2. Expansion of the Differential Operator. For any tensor field T of rank I 
on M, we have on Af~(p) the identity 

VgTi = -^D 43 T! - ^D 34 Tj + X ab D ab Tj 

= -D 43 Tj + ii? 43 [T] I + X ab D ab Tj 

= -f 4 {D 3 T) I + 2rf y Q Tj + ii? 43 [T]j + A ab Z? ab r 7 , 

where V 4 denotes the operator V L , where R 43 \T]i denotes the curvature quantity 
-D43T/ — D34T1, and where in the last step, we applied the identity D^L = 2rj a e a . 

By the definitions of the connections D, D, and t, we have the identities 
D ab Ti = D a {D b T) I - D^ a& T u 

fabTi = f a (f b T) i - fy ae Ti = D a {D b T) I - Dy^Tj, 
so by ©, we have the relation 

D ab T! = f ab T! - D^^^Tj = t ab T! - \xJiTt - ^Xab^Tj, 
As a result, we have the decomposition 



(30) 



□Tj = fa - y 4 (D 3 T)j + 2rft.Tr - - (tr*) tiTj 
-l(tr X )D 3 T I + ^R 43 [T} I . 



Thanks to Qian Wang and her comments for greatly simplifying the proof. 
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As a result, the quantity 

n ~ 

sf = E / ^ (m)/ * (m) / 

m=l ^ 

= E / A (m)I n g & m) ! + E / v4 (m,/ p (mc V/ J ^ $ 

where we have used |T]), can be expanded by ([50)1 as 



MSW 



n „ 

# = x: / a^ 



£$<"*> J - ^4 (Da*'" ), + 2^y Q $('"» / + -i? 43 I*'" 1 ']/ 



E ^ 



(m), 



- (trx) y 4 $ (m) / - 5 (trx) D 3 $ (m) , 



E 



^t^pW^/fl^Clj 



= a* + g£ + af + yf + a* + ^5 + a 6_ 

Moreover, we can expand 2lf as 

21" 






= £* + £ 2 e + € 3 C . 

4.3. Integrations by Parts. The main goal in principle is as follows. As of now, 
all the existing covariant derivatives thus far are acting on the $ (m) 's. We wish to 
move all covariant derivatives in directions tangent to Af~(p) (that is, L, e±, and e-i) 
away from the $ (m) 's. This will be accomplished through multiple applications of 
Lemmas [8] and [TO] We will then see that the transport equation (|T9]) will eliminate 
the worst remaining terms: those involving L-derivatives of the $ (m, 's. 
If we apply Corollary |H] to 21* , we obtain 



2te = E / -t A^^f^'j-f (logtf),4 <m)/ y a $ <m) / . 
Combining this with ([5]), then 

S£ + 3l? = y; [-/ fV-'f „*"»>/ + [ (rf-OA^t^j 

m=l L ^ • / ^ 

Applying Lemma [5] to 03 * yields 



m=l Jj ^> 



( m )J 



K (logi?) ■ (c a -?z a ) $<m) ' + (^ a Ca - T%) $ (m) , 



20 



ARICK SHAO 



n p 



m=l jAf < 
Likewise, we apply Lemma [TU] to 21^ to obtain 



21- 



E 

m— 1 



n 

El 



f 4 A (m>I D 3 & m) i 



A', 



(tT X )A^ 1 D 3 ^^ I 



A^Ds*™! 



A [m)I D 3 & 



m=l V°« J,S "o 

Taking note of the transport equation ([T9")) , then 






A (m)/ D 3 4 (m) / 



/ A w, J) 3 $ W f =£ e 1 +3 1 . 



Remark. TTie above computation also provides the justification for the transport 
equation (|19[) . Indeed, this is exactly what is needed to eliminate the terms involving 
integrals over M e of the £>3$ (m) 's. Recall from standard local energy estimates that 
L- derivatives on the null cone are generally the worst behaved. 

Remark. Recall also that [12] justified its transport equation as the condition 
needed to eliminate terms involving the distribution S'(u), where u is the "optical 
function" used throughout |12j . Note, however, that S'(u) is precisely —\ ■ D 3 S(u), 
which corresponds to an integral over the null cone of some terms, with the principal 
term being the L-derivative of the test function. 

We also apply Lemma fTUl to 21^: 



* = £ 



n 

E 

m— 1 
n 

E 

m— 1 

i4 i co5 



l -H 4 (trx) A^ 1 ^! + i (trx) lM <m)7 $ (m) J 



1 



(trx) (trx) A (m)/ $ (m) 7 



\ f (trx) ^'4'"'i - \ I (trx) A™'*'™', 



We can further expand 25^ using ©: 



Q3f + <8f 




y a R a + h\ + 4 ( tr x) (trx) - 2^ af> *a6 




+ Sl.G a ~-5 a -) i,w ' #M ' 




= Bj + b! 


+ < B 9 e +ml + m 2 e +m 3 f . 



A Ml&r, 



By the transport equation (|19l) . then 

®S + ®! = i E / ( tr *) p^/a""*™, = vl 

4 m,c=l- / - w ' e 
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If we apply Lemma [8] to <L\, we obtain 

n „ 



A/",: 



i ^ao 

m,c— 1 < 



£ / P<™> Q .//>>'V" 



q+^+S 



? + ^7 p - 



Moreover, applying Lemma fTOl to C^, we have 



771, C— 1 



£ e=-^ E / i Q ^ (C,J ^ <Cm, Q J / y4$ (m, / 

2 „_■, J AT, 



771, C=l e 

71 

+ E 



(tr X )P < "" ) 3jM (c)J $ (m) / +7 7 a P<™ ) aJ 7 ^ )J $< m, / 



ap(cm) IlW^H, 



AC V 2 



1-x 



y 4 P ( ™ ) 3,/A (c)J $ (m » 7 + ip (cm) 3J / y4^ <C,J * (m> / 



E (If p lm WA^ j i^ I -^f p^w 

m,c=l \ ^ S « ^ 5 "0 



A< C > J $('") / 



= £° + q + qjf + €' e + Z\ + 3 Z . 
Applying the transport equation (fl"9)) . we have 



n 1 - 

£" + <= E 1/ (tr^P^a/A'""*™, 
m ,c=l 4 •'AT. 

n , . 

+ E 7 / P (cd) 4 J K P (d '" , 3K / A (c)J $ <m) / 



m,c,d— 1 
}4 i m5 



As of now, we have the rather unsightly expansion 

(3i) sf = #7 + ^7 P + g< + yf + <&l + Bf + ® J + Q3^ + < + erf + an, 1 + gj£ 

+ Ml + & + «# + <p3 + qj4 + ^5 + £ + £l + £ 2 + 3 o + jl + 3 2 

4.4. The Error Terms. The next step is to aggregate the terms in (|3i"j) into error 
terms corresponding to those comprising (p?Tj) . First, recalling ([5]), we see that 

m=l ^ 
n « 

Recalling the mass aspect function /i defined in ((7j) , we have 

n „ 

m\ + an* + ml + m 4 e + wi 5 e = Y, / M Cm)/ $ <m) / = 5?- 

1 JAC 
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Next, by ©, 
hence we obtain 

n „ 

< + <+ ^ + <£ 2 +<$+<#+ ^ = ^ / v l«*/A™ J & m > I =ft. 

m , c = l We 

Finally, we obtain the more manageable equation 

(32) yf = s £ v + & e + $J P + $» + r e + yf + £° + £* + £ 2 + 3° + 3 1 + a 2 . 

Note that: 

• The "#" -terms are integrals over the regular past null cone segment A/" e . 

• The "£"-terms are the vertex limit terms, expressed as integrals over S e . 

• The "3" -terms are the initial data terms, expressed as integrals over S Vo . 
The limits as e \ of the "3"" -terms in (|32p are clear - these involve simply replacing 
the integral over J\f e by the same integral over J\f~{p;vo). These account for the 
terms $(p;v ), ("^(pjuo), and £ 2 (p;vo) in TheoremlTl Moreover, the "3"-terms in 
([52")) correspond exactly to 3(p;«o) in TheoremlT] 

As a result, the limit as e \ of (f3"2"| evaluates to 

(33) - lim (£° + £i + £ 2 ) = J (p; u ) + (S 1 (p; u ) + (£ 2 (p; i; ) + 3 (p; t, ) . 

Now, it remains only to compute the limits on the left-hand side of Q33p. 

4.5. Vertex Limits. To compute the limits of the "£"-terms, we first obtain the 
initial values of the integrands using Proposition |6l and then we apply the integral 
limit property (|16p to each term. 

We begin with £ £ , which we write as 

n „ n „ 

m=l^ 5 ' m=l ''• S « 

Along each null generator 7^, y G S 2 , we have the limits 

As a result, by P^|) . 

lim£* =0. 

Similarly, for £ 2 , we have 

1 " r 

lim£ 2 = - V lime- 1 / pM ^W^W = . 

m,c— 1 ° e 

Next, for £ £ , we write 

1 " r 

£ ME^ 2 / [ e (trx)]S (m,/ <J> ( ™' 7 . 
m=l ^ 5 < 

Equation (TT5|) implies e(trx) converges to — 2i9q 1 at the vertex, hence by (|16|) . 

* n „ n 

lim£° = 2^0 E / ("^o" 1 ) ^ $ (m, /| P ^ = -4^o £ J(m>/ *'"%• 

C m=lV s2 m=l 
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Therefore, we obtain the desired limit 

n 

- lim (£° + Z\ + Zl) = 4^o £ J(m>/ $(m ^l P ■ 

m— 1 

As a result, equation ([55]) becomes (j^Tj) . and the proof of ([2"Tj) - ([2"5j) is complete. 
Lastly, for the alternate representation (1261) of the main error terms, we simply 
apply Corollary [9] to (l23l) . This concludes the proof of Theorem [71 

Remark. We can also apply the method of proof in this paper to derive the repre- 
sentation formula KR of [12] - The calculation is then significantly simplified, since 
the term 2tf in the proof now vanishes, as do the terms resulting from 21® . 

5. Generalization to Vector Bundles 

The representation formula KR was further extended in [12l Thm. 4.1] to the 
case of covariant wave equations on sections of arbitrary vector bundles. For con- 
venience, we denote this extended parametrix by KRV. We can produce a similar 
extension of Theorem [7] to sections of vector bundles, using a process analogous 
to that of [12j . We summarize this method in this section, and we connect this 
extended version to Theorem [JJ as well as to KRV. The final result represents a 
culmination of the ideas which led to KR, KRV, and Theorem [JJ 

5.1. General Ideas. The extension from KR to KRV is essential for the main 
application presented in (T2] : a gauge-invariant proof of the L°°-estimates for the 
Yang-Mills curvature required in Eardley and Moncrief 's classical global existence 
result for the Yang-Mills equations in [5]. o Recall that [5] applied the standard 
representation formula for the scalar wave equation in Minkowski space and, as 
a result, required the Cronstrom gauge condition in order to control the gauge 
potential. The same was also true for its curved spacetime analogue [3j. Using 
KRV, however, one can avoid the gauge potential altogether and hence bypass the 
need for a favorable gauge condition. 

Note that the parametrix KR itself is not directly applicable to the Yang-Mills 
setting. In order to achieve this gauge-invariance property, one requires a version 
of the representation formula suitable for Lie algebra- valued tensor fields and the 
associated gauge covariant derivatives. Not only does the extended parametrix 
KRV of [12] Thm. 4.1] address this issue, but its derivation is also completely 
straightforward in light of the proof of KR. 

The key observations behind this extension concern the natures of the operations 
used in the derivation of KR. These can be summarized by the following statement: 

• On the past null cone Af~(p) for which we derive the parametrix, the 
(mixed) covariant derivatives on A/" - (p) satisfy natural Leibniz rules and 
are compatible with both the spacetime metric g and the induced horizon- 
tal metrics A on M~{p). In particular, this legitimizes the integration by 
parts operations fundamental to the derivation of KR. 
In other words, if covariant derivatives on a vector bundle behave in a manner 
analogous to the above statement, then the proof of KR can be directly adapted 
to prove KRV. In particular, the gauge covariant derivatives on Lie algebra-valued 
tensor fields in the Yang-Mills setting satisfy a version of this property, as long 



•*In fact, the papers [UE dealt with the more general Yang-Mills-Higgs equations. 
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as the associated metrics are modified using an appropriate positive-definite scalar 
product on the Lie algebra under consideration. 

It is unsurprising, then, that an analogous extension can be made to the repre- 
sentation formula of TheoremUJ The key observations behind this extension remain 
the same as before, so once the framework is set, the extension becomes completely 
straightforward. Next, we will briefly summarize the basic constructions for this 
extension of Theorem [7j 

5.2. Basic Constructions. Let V denote a vector bundle over M, and let V and 
(•, •) denote a connection and a bundle metric on V, respectively. We additionally 
stipulate that T> and (•,•) satisfy the following compatibility condition: for any 
sections S, T G FV, we have the Leibniz-type identity 

D a (S, T) = (V a S, T) + (S, V a T) . 

Recall that TV denotes the space of all smooth sections of V. 

Next, we let V denote the restriction of V to 7V~ (p), and we note that T> induces 
a connection I? on V. We will also require the following constructions: 

• We can define various "mixed bundles" over Af~ (p) , expressed as 

T k T r V = T k Af- (p) ® T r Af~ (p) ® V, k, r > 0, 

where the right-hand side denotes tensor products of vector bundles. 

• Moreover, we can naturally construct "mixed covariant derivatives" V on 
the above mixed bundles, which behave like V on the "horizontal" com- 
ponents, D on the "extrinsic" components, and V on V. In particular, for 
a decomposable element A eg) B ® C € TT_ T V and a vector field X on 
7V~ (p) , we define the mixed covariant derivative by 

f x {A ® B (g) C) = f x A ig) B <g> C + A ® D X B <g> C + A <8> B ® V X C. 

• The metrics (•, •), g, and 7 induce pairings for the mixed bundles, which we 
also denote by (•, •), and which, similar to the original bundle metric, satisfy 
compatibility conditions with respect to the above mixed connections. 

• We can define the V-curvature "R. to describe commutations of two 2?- 
derivatives. More specifically, we can define for any T G TV and vector 
fields X, Y on M the quantity 1Z X y [T] S TV by the formula 

TZ XY [T] = V x (V Y T) - V Y (V X T) - T> [x ,y)T. 

We can recover the setting of KR by taking V to be the bundle T r M of all 
tensors on M of total rank r, with T> = D the connection on T r M induced by the 
Levi-Civita connection and (•, •) = g the full metric contraction operation on T r M. 
Note that V and (•, •), as given here, are indeed compatible. For this special case, 
the above mixed bundles and their associated connections and metrics coincide with 
the corresponding "mixed" objects defined in previous sections. This construction 
of mixed bundles and derivatives along with the compatibility of the associated 
connections and metrics are the driving forces behind the extension of Theorem [7] 

We will denote the object satisfying the covariant wave equation on V to be 
$ G FV, and we denote the associated nonlinearity by \& G TV. It remains now 
to describe the first-order coefficients of the wave equation, which we denote by P. 
In this setting, we can define P to be a "(V ® V)-valued 1-form", that is, a section 
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of the vector bundle T*M V V, where T*M is the cotangent bundle of M. In 
particular, the restriction of any component P a to 7V~(p) is a section of V V. 
Using the bundle metric (•, •), we can define the natural bilinear pairings 

I-,-) : r(V0V) xrv^rv, (-,-| : fvxf(V0V) -»rv, 

given for decomposable elements by 

\A®B,C) = A-(B,C), (A,B®C\ = (A,B)-C, A,B,CeTV, 
and extended accordingly to general elements. We will also need the maps 

(• | • | •) : rv x r (v v) x rv -»• c°° (Af- ( P )) , 
{•,•}: r (v v) x r (v v) -»• r (v v) , 

defined for decomposable elements by 

(A I B C I D) = (A, B) (C, D), {A® B,C ® D} = (B,C) ■ A® D, 

and extended as multilinear functions to general elements. 

For example, consider the standard case V = T r M of KR. Let S,T € TV, let 
A,B^ r(V V), and let /, J, K denote collections of r extrinsic indices. Then, the 
above operations are given explicitly via tensorial contractions, i.e., 

\A,T) I =A I J T J , (T,A\ J = T I A I J , 
(S \A\T) = SjA ij Tj, {A, B} u = A IK B K , 7 . 

We have a natural connection I? on V V induced from the connection T> on V, 
given by the Leibniz-type formula 

V (A B) = VA B + A VB 

for the decomposable elements and extended linearly to general sections. Then, we 
can once again define mixed bundles and mixed covariant derivatives ^ correspond- 
ing to the bundles V V. It is a simple exercise to determine that the expected 
compatibility conditions hold for |-, •), (•, -|, (• | • •), and {•, •} with respect to the 
various connections and covariant derivatives. We omit the rather tedious details. 

5.3. The Generalized Equations. The covariant wave equation under consider- 
ation is now given by 

(34) D^$ + \P, T&) = gTV^* + \P„, V^) = *, 

where $, $ e TV and P e F(V V). The precise definition of \Jf can be a bit 
subtle. Here, the inner covariant derivative T>p is that of the connection on FV, 
but the outer covariant derivative V a is the induced connection on T(T*M 0TV), 
which acts like D on T*M and like V on V. In other words, for vector fields X, Y 
on M, we define the "second derivative" T>xy$ by 



v XY <§> = v x 


(Vy$) - V DxY <P. 




Let B € TV satisfy the system of transport equations 




5) t f B = - 1 - 


0(trx)-| 


B + ^(B,P 4 \, 


B l 



J. 
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where J is element of the fiber at p of V. We also define A = f~ l B G TV, as well 
as the coefficients v £ T(V £S> V), given by 

V = ^Pa + ^P? + CPa + \ (tl x) P 4 + J (tr *) P 3 + ^ {P 4 , P 3 } , 

We have indexed with respect to arbitrary /-adapted null frames on Af~(p). 
We can now finally state our extension of Theorem [7j 

Theorem 11. Suppose V is a vector bundle over M , with compatible connection 
and bundle metric T> and (•,•). We assume all the objects derived from V, T>, and 
(•,•) as defined in this section, and we define <&, $!, P, B, A, v as in the above 
development. If p £ M, and if the past regular null cone J\f~(p) satisfy the same 
assumptions as in Theorem^ then 

(36) 47n?o (J, $| p ) = 3 {p; v ) + (B 1 (p; v ) + <E 2 (p; v ) + 3 (p; wo) , 

ff(p;«o) = -/ (A*), 

JN'-fD-.Vn) 



where 

$(p;vo) = - 

I 

>N-(p;v ) 



e 1 (p; v ) = - [ \(H a A, V a $) + (c a - ? 7 a ) (faA $ 

JH-(p;v ) LX ' ^ 

£ 2 (p; «o) = / m (A, *) + ~ / (A, n 43 [*]) 

JAf-(p;v ) Z JAr-(p;v ) 



Af-(p;v ) 



(ll/ a A | P„ | *) + <A | i/ | *>^ 
3(p;v ) = -\[ (trx)(A,*)- / <A,X> 3 *>-~/ (A | P 3 I *) ■ 

*» Olln «* Oil,-, «* OlJfi 



/s. 

Again, the error terms <£}-(p;vo) can be alternately expressed as 

€ 1 (p;vo)= [ (jA + 2Cf a A,<p) 

JM-(p;v a ) X ' 

We can derive Theorem [TT] in a manner completely analogous to Theorem [7J We 
simply replace the compatibility properties of Z? and y in the tensorial case by the 
abstract compatibility properties described above. 

5.4. Reduction to Previous Cases. Note that we can obtain KRV trivially 
from Theorem 1111 simply by taking P = 0. Next, we will show how to Theorem [7] 
can be recovered from Theorem 111! 

For this, we consider the direct sum vector bundle 

V = T r<1) M©...©T r( " ) M, 
with the associated metric and connection on V given naturally by 

n 

(37) ((SW. . .,S<»>) , (TW . . . ,T<">)> = £ 5 (S<™>,T<»>) , 

771—1 

P(S (1) ,...,S ( ' l) ) = (I?S (1 \..., J DS ( " , )■ 
We then define the fundamental objects 

$ = ($ (1) ,..., $<")), tf =(tf (1) ,...,* (B) ), 
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and we similarly define P from the component objects _P< mc >. With these specific 
definitions, then the covariant wave equation (f34|) coincides with ([TJ), and formula 
of Theorem QTJ reduces to Theorem [7] 

Note that the vector bundle formalism of (|36|) automatically encapsulates the 
fact that Theorem [7] deals with a system of tensor wave equations. As a result, 
the parametrix KRV can also deal with systems of covariant wave equations in 
this fashion. What is entirely new with respect to 12 is the presence of first-order 
terms, the resulting coupling, and the weakening of the assumptions. 

Remark. Alternatively, one could possibly reinterpret the first-order coefficients P 
by modifying the given connection T> to also take into account contributions from P. 
This would bring us closer to the case o/KRV. However, such modified connections 
would generally fail to be compatible with the given metric, hence one would need 
to deal with error terms resulting from this. 

5.5. Applications to General Relativity: Einstein- Yang-Mills Spacetimes. 

We conclude this paper by briefly summarizing a couple of applications for which 
Theorem ITTI could be helpful. First, we can apply Theorem ITTI to treat the analogue 
of Theorem [5] for the Einstein- Yang-Mills setting. This is quite analogous to the 
Einstein-Maxwell problem discussed in [16] , except that the matter field "$ is now a 
0- valued 2-form, where g is the Lie algebra associated with the Yang-Mills theory. 
In other words, we can think of $ as a section of the bundle T 2 M <g) g. 

Recall also that we can define "gauge covariant" derivatives on such g-valued 
tensor bundles, which we denote by T>. Moreover, if we assume as is standard that 
q admits an invariant positive-definite scalar product (■,•), then we can combine 
(-,-) and g to define bundle metrics on such g- valued tensor bundles which are 
compatible with T>. Hence, we retrieve a setting that can be described using the 
language of Theorem ITTI 

Similar to the Einstein-Maxwell setting, we again have a system of two covariant 
wave equations in the Einstein- Yang-Mills analogue, one for the spacetime curvature 
R, and one for the derivative of the Yang-Mills curvature, T>$ £ T(T 3 M <g) g). 
Algebraically, these equations are analogous in form to the system (|28|) treated in 
the Einstein-Maxwell setting. Unlike in [16] , we cannot apply Theorem [7] here, 
since the wave equation satisfied by T>$ involve gauge covariant derivatives rather 
than the standard tensorial derivatives, lj We can, however, describe this setting 
in terms of Theorem I 111 In fact, we simply take 

V = T 4 M © (T 2 M ® fl) , 

and we construct a natural connection and bundle metric on V from the above 
developments using the same method as in (I37|) . Consequently, the system of wave 
equations described above is now of the form (|34|) , and we can apply Theorem [TTJ 
in order to derive the necessary uniform bounds for R and T>$. 

5.6. Applications to General Relativity: Time Foliations. In [20j [21], and 

to a lesser extent in [16] , one requires the observation that the second fundamental 
form k of the given time foliation satisfies a covariant wave equation. This can be 
cleanly described in terms of the language of this section. Furthermore, with this 



In terms of the usual tensorial covariant derivative, such a wave equation for 5 would of 
course fail to be gauge-independent. 
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characterization, one can more easily derive the energy estimates for k associated 
with this wave equation, which are crucial in the above works. 

In this setting, the spacetime M is given as a smooth 1-parameter family of 
spacelike hypersurfaces {S r }, and k denotes the second fundamental forms on the 
E r 's. Then, the bundle V under consideration here is the "horizontal bundle" of 
all symmetric 2-tensors in M which are tangent to the S T 's; note that k £ TV. For 
the bundle metric, we can simply consider the pullback 7 of g to the S T 's, i.e., the 
induced metrics of the S T 's. The compatible connection on V can be defined in the 
same manner as for the connections on the horizontal bundles of past null cones. 

From the Einstein equations, one can see that k satisfies a covariant wave equa- 
tion of the form (|34[) . As a result, we can apply Theorem ITTI to k in this setting. 
In fact, the application of a representation formula to k is an essential component 
in the arguments of [20l [21] . Note that in contrast to the previous example and to 
Theorem [3 the bundle metric in this case is in fact positive-definite. 

Recall that in the case of scalar wave equations, one can derive energy estimates 
using energy-momentum tensors in a completely standard fashion. Moreover, this 
methodology extends directly to the setting of an arbitrary vector bundle along 
with a connection and a compatible positive-definite bundle metric. Therefore, by 
expressing k in the form (|34[) , then one can immediately derive associated global 
and local energy estimates for k. 
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